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Abstract. - Using the Chalker-Coddington network model as a drastically simplified, but uni- 
versal model of integer quantum Hall physics, we investigate the plateau-to-insulator transition 
at strong magnetic field by means of a real-space renormalization approach. Our results suggest 
that for a fully quantum coherent situation, the quantized Hall insulator with i?H ~ h/e 2 is 
observed up to _Rl ~ 25/i/e 2 when studying the most probable value of the distribution func- 
tion P(Rh)- Upon further increasing J?l — * oo the Hall insulator with diverging Hall resistance 
Rh (X J?l i s seen. The crossover between these two regimes depends on the precise nature of 
the averaging procedure. 



Introduction. - The integer quantum Hall (QH) transitions are described well in terms 
of a series of delocalization-localization transitions of the electron wavefunction [1]. These 
universal plateau-plateau transitions are accompanied by a power-law divergence e _,y of the 
electronic localization length £, where e defines the distance to the transition for a suitable 
controlled parameter, e.g. the electron energy [2]. Similarly, it is now conclusively established 
that plateau-plateau and insulator-plateau transitions exhibit the same critical behavior [3-6] . 

However, the value of the Hall resistance Rh in this insulating phase (at large magnetic 
field) is still rather controversial. Various experiments have found that _Rh remains very close 
to its quantized value h/e 2 even deep in the insulating regime [4-6] with longitudinal resistance 
i?L > h/e 2 . This scenario has been dubbed the quantized Hall insulator. On the other hand, 
theoretical predictions show that a diverging i?H should be expected, i.e., Rh cx i?£ [7,8]. 
This Hall insulator is to be expected at strong disorder or strong magnetic fields. 

In fully quantum coherent transport measurements such as in mesoscopic devices at low 
temperature, the results clearly show the paramount influence of quantum interference and 
the measured quantities fluctuate strongly [9]. At magnetic field B — 0, the universal con- 
ductance fluctuations provide the most famous example [10]. For the QH situation, similarly 
reproducible and pronounced fluctuations have been observed previously [11-14], although no 
complete understanding of their behavior has yet emerged. Thus for quantum coherent QH 
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Fig. 1 - RG unit of 5 SP's (full circles) and equipotential lines (arrows) on a square lattice with 
Eq. £Q| combining them into a super-SP. $i, . . . , $4 are the phases acquired by an electron drifting 
along the contours. Arrows show the electron motion for one B direction and should be reversed for 
—B. The small solid dots indicate the position of the weakly-coupled voltage probes. The thin lines 
denote the different possibilities of measuring voltage drops in the present structure for both B-field 
directions. 

samples, the average values J2l and become meaningless unless their full distributions are 
being taken into account [15]. Consequently, different averaging procedures may then yield 
quite different estimates for J?l and -Ru- 
in the present manuscript, we revisit the Hall insulating regime by using a real-space renor- 
malization group (RG) approach to the Chalker-Coddington (CC) network model, recently 
used successfully to study the energy-level statistics at the QH transition as well as the influ- 
ence of long-range correlations in the disorder potential close to the transition [16]. We will 
show when considering the experimentally relevant most probable value as estimate for i?n, 
we find quantized i?n = h/e 2 for i?L < 10/i/e 2 and increasing values for larger i?L thereby 
reconciling the experimental results of a quantized i?n available up to i?L = 8/1/e 2 [5,6,14] 
with the theoretical predictions of the Hall insulator. 

Model and RG approach to the QH situation. - The CC model is based on the microscopic 
picture of electron motion in a strong magnetic field and a smooth disorder potential [17] when 
only the semiclassical trajectories of the guiding center of the cyclotron orbit are important. 
Assigning these trajectories to links and considering saddle points (SP's) at which different 
trajectories come closer than the Larmor radius as nodes a chiral network can be constructed. 

We now apply a real-space RG approach [11,18] to the CC network [8,12,16,19]. The 
RG unit we use is extracted from a CC network on a regular 2D square lattice as shown in 
Fig. ^ A super-SP consists of five original SP's by analogy to the RG unit employed for 
the 2D bond percolation problem [20]. Each SP can be described by two matrix equations 
relating the wavefunction amplitudes in incoming and outgoing channels. Between the SP's 
an electron travels along equipotential lines, and accumulates a certain Aharonov-Bohm phase 
4>. Different phases are uncorrelated, which reflects the randomness of the original potential 
landscape. This results in a system of ten linear equations, the solution of which yields the 
expression for the transmission coefficient of the super-SP [11] 

t , = tit 5 (r 2 r 3 r 4 e t * 3 - 1) + ^t 4 e'^ 3 +^)(r 1 r 3 r 5 e-^ 1 - 1) + t 3 (t 2 t 5 e , *» + hhe^) 

(ra - r 2 r 4 e»*>)(r 3 - nr 5 e^) + (i 3 - t 4 t 5 e***)(t 3 - ht 2 e^) ' U 

Here and r, = (1 — if) 1 / 2 are, respectively, the transmission and reflection coefficients of 
the constituting SP's, <&j are the phases accumulated along the closed loops (see Fig. HJ. 



Philipp Cain and Rudolf A. Romer: Fluctuating Hall resistance defeats the quantized Hall insulator3| 



Equation Q is the RG transformation, which allows one to generate (after averaging over 
$j) the distribution P(t') of the transmission coefficients of super-SP's using the distribution 
P(t) of the transmission coefficients of the original SP's. Since the transmission coefficients 
of the original SP's depend on the electron energy e, the fact that derealization occurs at 
e = implies that a certain distribution, P c (t) — with P c {t 2 ) being symmetric with respect 
to t 2 = \ — is a fixed point (FP) of the RG transformation The distribution P C (G) of 
the dimcnsionless two-point conductance G can be obtained from the relation G = t 2 , so that 
P C {G) = P c (t)/2t. We remark that the classical percolation limit [20] is obtained when the 
classical values t € {0, 1} are used in the RG. 

Let us now make contact with the resistances. The dimensionless longitudinal resistance 
i?L can be computed from G via 

\_rf_ _ 1 - \t\ 2 _ 1_ 
1*F ~~ \tP ~~ G 



R L - TTT2 ~ ui2 = 7* ~ 1 = ifet ~ 1 (2) 



with the dimcnsionless 2-terminal resistance i?2t ■ In order to study the Hall resistance for the 
four-terminal CC model, we study the difference in (normalized) chemical potentials ^15 and 
/Z42 at the positions indicated in Fig. Q] [7, 8]. In order to remove any cross-current amplitudes, 
the voltage is computed, analogously to the experimental situation, by subtracting opposite 
magnetic field directions. The Hall voltage is then defined as 

Un=~ {[Mi5(5) - fUa{B)] - \pl 16 {-B) - M4a(-B)]} • (3) 
Consequently, the dimensionless Hall resistance is given by 

*H = ^ • (4) 

Numerical RG procedure. - In order to find the FP distributions P c (Rl) an d P c (Ru) at 
the QH transition, we start from an appropriate initial distribution of transmission coefficients, 
Po(t). The RG procedure results in a broad distribution P C (G) which is peaked at G > and 
G < 1 [16]. Any RG flow away from this instable FP towards the insulating regimes will 
result in a further increase of one of these peaks. Therefore, in order to reliably model the 
insulating regime, an accurate procedure for reproducing these peaks is needed. In [16], we 
had also shown that a distribution Q(z) of dimensionless SP heights is related to P(G) via 
Q(z) = P{G){dG/dz) = \ cosh" 2 (z/2)P [{e z + l)" 1 ] and t = {e z + I)- 1 ' 2 . Thus the peaks in 
P(G) transform into long but quickly decaying tails in Q(z) and it is numerically much better 
to perform the RG for SP heights z when the physics away from the FP is studied. 

We discretize the distribution Q(z) in at least 6000 bins such that the bin width is typically 
0.01. Since z e] — 00, oo[, we have to include lower and upper cut-off SP heights such that 
z £ [ziowjZup]. We use zi ow — —20, z up = 40 for perturbations towards positive z and vice 
versa for negative perturbations. This corresponds to transmission amplitudes t(z\ ow ) w 10~ 5 , 
t(z up ) wl - 10~ 9 for positive perturbations and vice versa for negative perturbations. From 
Qo(z), we obtain Zj, i = 1, . . . , 5, compute the associated t, and substitute them into the RG 
transformation 0J. The phases j = 1, . . . , 4 are chosen randomly from the interval $j £ 
[0, 27t]. In this way we calculate at least 10 8 super-transmission coefficients t' and associated z'. 
At the next step we repeat the procedure using Qi as an initial distribution. We assume that 
the iteration process has converged when the mean square deviation J dt[Q n (t) — Q n -i(t)] 2 
of the distribution Q n and its predecessor Q n -\ deviate by less than 10~ 4 . We check that 
the values of zi ow ,u P do not influence our results and that all the previous results at the 
QH transition as in Refs. [16] are reproduced. The full width at half maximum of the FP 
distribution Q c (z) is about 5 [16]. 
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FP distributions. In Figs. [3 and [3 we show the FP distributions for i?L and i?n, 
respectively. As seen previously [7], both distributions are manifestly non-Gaussian with 
rather long tails. Indeed, P(Rl) can be fitted by a log-normal distribution as indicated in 
Fig- 13 [7]- The FP distribution P(Rn) exhibits a pronounced peak at Ru = 1 (h/e 2 in SI 
units) and a small shoulder around i?H = 0. The small weight for negative Hall resistances 
shows that in the coherent RG structure chosen, the large fluctuations in P(t) can result in 
an apparent reversal of the Hall voltage. 

The plateau-insulator transition. - In order to model the transition into the insulating 
regime, we now shift the initial distribution Qo(z) — > Qo(z — zq) by a small zq. For zq < and 
zq > 0, the RG flow will then drive the distributions into the insulating, G — ► 0, and plateau, 
G — ► 1, regimes, respectively. In Figs.|3and|3 we show the resulting distributions after many 
RG steps. For P(Rl), we see in Fig. [3 that the flow towards G — > 1 results in a decrease 
of large i?L events, whereas conversely, the regime G — > leads to an increase in large i?L 
values and a decrease of the maximum value in P(R^). For P(i?n), Fig. [3 shows that the 
plateau regime G — ► 1 gives a highly singular peak at the dimensionless quantized Hall value 
1, corresponding to a perfect Hall plateau. On the other hand, the insulating regime G — > 
shows an increase of weight in the tails of P(Rjj) and the eventual obliteration of any central 
peak. 

In order to extract an averaged i?L and i?n from these non-standard distributions, we 
should now select an appropriate mean (•) that characterizes and captures the essential physics 
and allows comparison with the experimental data. The precise operational definition is 
also important as it corresponds to different possible experimental setups. Therefore, we 
consider several means: (i) arithmetic (i?) ar i = ^2iRi/N, (ii) geometric/typical (i?)t yp = 
expj^i lo.Ri/N, (iii) median (central value) (i?) mo d, where N denotes the number of samples 
(> 10 8 ) in each case. The median and the typical mean (and their variances [8]) are less 



Fig. 2 - Distribution of the longitudinal resistance _Rl- The □ symbols indicate the FP distribution, 
the solid line is a fit to a log-normal distribution and the dot-dashed and dashed lines show the 
distribution after n = 15 RG iterations into the conductance regimes G(zo < 0) — > and G(zo > 
0) — > 1. Only every 5th data point is shown for P(Rl). 

Fig. 3 - Distribution of the Hall resistance Ru. The FP distribution has been shaded to zero. For 
the G(zo < 0) — » and G(zo > 0) — > 1 regimes, we have indicated the distributions after n — 16 RG 
iterations by bold dashed and dot-dashed lines. 
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sensitive to extreme values than other means (such as, e.g. root-mean-square and harmonic 
mean) and this makes them a better measure for highly skewed and long-tailed distributions 
such as P(.Rl) and P(i?n) in the insulating regime. In the plateau regimes, the distributions 
are less skewed, particularly for i?n and the difference in the means becomes less important. 

We are left with determining in which operational order to apply the averaging procedure. 
For i?L, as measured via Ohm's law as a ratio, it is obvious that the appropriate average 
should be i?L = 7^y — 1 (and not (^} — I). For Rn, the situation is less straightforward due 
to the definition of Un in ©. A simple average is (£/h), i.e. using the appropriate P(i?n)- 
Similar to the experimental procedure, we can also estimate C/h via (^15-/^42) for each B field 
direction separately. In Ref. [8], it has been suggested that a more appropriate average (C/h)* 
can be constructed from (/U15} — (^42)- This later procedure corresponds to measuring the 
voltage drop between positions fii$ and /i42 in Fig. ^ by separately measuring the individual 
voltages with respect to ground and then recombining them. 
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Fig. 4 - Dependence of averaged Rh on averaged Rl for various means. The dashed line describes 
the divergence Rh oc 7?£ of the typical mean {i?H)typ, the median (7?H)m e d an d the most probable 
value of P(Rh)- The variance of the arithmetic mean diverges as _Rl — * oo and (-RH)ari is no longer 
a useful characteristic of P(Ru)- The horizontal dotted line indicates _Rh = 1. 

Fig. 5 - Plot of the same data as in Fig. 2] but on a linear scale and for experimentally accessible 
resistance values. Here the almost quantized behavior for the most probable value of P(Rh) becomes 
even more pronounced. The horizontal dotted line indicates _Rh = 1. 



In Figs.Q]and[S]we show the resulting dependence Rh(Rl) when these averaging definitions 
are being used. These figures are our main result. In the plateau regime i?L < 1, the nearly 
constant behavior of all i?n averages is as expected. For i?L — > oo, we get a divergent Ru 
when using the median and typical means as suggested by Refs. [7,8] to compute both i?L and 
i?H- This divergence can be captured by a power-law i?H oc i?£ with k « 0.26. The arithmetic 
mean for large i?L 3> 1 quickly becomes instable and no useful information can be inferred. 

Reducing the information to the experimentally more relevant resistance regime of a few 
times h/e 2 , we replot the i?L and i?H data in Fig. In addition to the three means above, 
we also show the behavior of the most-probable value at which P(i?n) has a maximum. 
This estimate Ph(Pl) appears relevant in the experimental setup where I0 8 different samples 
cannot be easily measured and the full distribution functions cannot be constructed in similar 
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Fig. 6 - Dependence of the localization lengths £r, and £h on the initial shift zo using the ansatz (8) 
of [7]. We obtain a reasonable estimate of the critical exponent v. 

Fig. 7 - The semi circle law for different means as in Fig. The most probable values (o) show 
perfect semi-circle (solid line) behaviour. The thin dotted line denotes <tl = 1/2. 



detail. Most importantly, for the range of i?L values shown in Fig- El the value of i?H deviates 
only slightly from its quantized value 1 at the transition and Rl + i?H ~ R2t [14]. Therefore, 
the experimental estimate of i?H appears to support the notion of the quantized Hall insulator. 
Indeed, the deviations from 1 are less than 10% until i?L ~ 25. However, going back to Fig. 
01 we see that in the strongly insulating regime Rl — ► 00 also Rn diverges with a power-law 
that is well-described by (R)* yp . We emphasize that fluctuations in i?u for large i?L ^> 10 5 
are due to numerical inaccuracies in P(i?u) arid decrease upon further increasing the number 
of samples. 

The results for i?L and in the localized regime can be very well described by an 
exponential scaling function with finite-size correction [7] R h< n(2 n ,z) oc 2? n exp[2 n ^ H (z)]. 

Plotting £l,h as a function of small perturbation zo, we find £l,h(zo) oc z lyL,H with j/l.h ~ 2.35 
as shown in Fig. Thus we recover the universal divergence of the localization length £ 
even when using the RG in the insulating regime [22]. An equally reliable estimate of the 
irrelevant exponent 7 appears not possible for our data. As usual, the 4-terminal resistances 
can be converted into the respective conductances via ctl,h = Rl.h/ {Rj, + i?n)- For these 
conductances, one expects the semi-circle law o\ + (or — 1/2) 2 = (1/2) 2 to hold [23]. In Fig. 
{7\ we show that the most-probable values capture the overall shape and symmetry properties 
best [5, 14], the other averages show pronounced deviations. We also note that the relation 
^l(zo) = 1/Rl(— zo) is obeyed by all means [21]. This is a consequence of the reflection 
symmetry of P(G) [16]. 

Conclusion. - We have shown in a quantum coherent calculation that the insulating Hall 
behavior of i?n(-RL — * 00) is dominated by the power-law divergence i?H oc R£. However, up 
to previously experimentally reached i?L < 10 [5,6,14], the deviation from a quantized i?u = 1 
is very small and the onset of the divergence for i?L 3> 10 is yet to be explored. 
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